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ABSTRACT 

Phase referencing is a standard calibration procedure in radio interferometry. It allows to detect weak sources by using quasi- 
simultaneous observations of closeby sources acting as calibrators. Therefore, it is assumed that, for each antenna, the optical paths 
of the signals from both sources are similar. However, atmospheric turbulence may introduce strong differences in the optical paths of 
the signals and affect, or even waste, phase referencing for cases of relatively large calibrator-to-target separations and/or bad weather. 
The situation is similar in wide-field observations, since the random deformations of the images, mostly caused by atmospheric tur- 
bulence, have essentially the same origin as the random astrometric variations of phase-referenced sources with respect to the phase 
center of their calibrators. In this paper, we present the results of a Monte Carlo study of the astrometric precision and sensitivity of 
an interferometric array (a realization of the Square Kilometre Array, SKA) in phase-referenced and wide-field observations. These 
simulations can be extrapolated to other arrays by applying the corresponding corrections. We consider several effects from the turbu- 
lent atmosphere (i.e., ionosphere and wet component of the troposphere) and also from the antenna receivers. We study the changes in 
dynamic range and astrometric precision as a function of observing frequency, source separation, and strength of the turbulence. We 
find that, for frequencies between 1 and 10 GHz, it is possible to obtain images with high fidelity, although the atmosphere strongly 
limits the sensitivity of the instrument compared to the case with no atmosphere. Outside this frequency window, the dynamic range 
of the images and the accuracy of the source positions decrease. We also find that, even if a good model of the atmospheric turbulence 
(with an accuracy of 99%) is used in the imaging, residual effects from the turbulence can still limit the dynamic ranges of deep, 
high-contrast (10 5 - 10 6 ), images. 

Key words, atmospheric effects - techniques: high angular resolution - techniques: interferometric - telescopes: SKA 



1. Introduction 

It is well-known that ground-based astronomical observations 
are affected by the atmosphere. Changes in the atmospheric 
opacity produce a bias in the source flux density, while changes 
in the refraction index distort the shape of the electromagnetic 
frontwave of the source. Such a distortion translates into a de- 
formation of the observed source structure and/or a variation of 
the relative positions of all sources observed in a given field. In 
the case of astronomical devices based on interferometry, atmo- 
spheric effects can be well modelled if the atmosphere above 
each element of the interferometer (hereafter, station) remains 
unchanged over the whole portion of the sky being observed. 
In such cases, the observed visibilities can be calibrated using 
station-based algorithms, which are relatively simple and com- 
putationally inexpensive (e.g. Readhead & Wilkinson 1978). 

However, when the spatial variations of the atmosphere are 
significant within the observed portion of the sky, as it happens 
if there is atmospheric turbulence, the opacity and dispersive 
effects cannot be modeled as a single time-dependent station- 
based complex gain over the field of view. Unless more com- 
plicated calibration algorithms are used (e.g., van der Tol, Jeffs, 
& van der Veen 120071 1. the effect of these errors on the image 
are difficult to correct. In this paper, we report on a study of the 
effects that a turbulent atmosphere may introduce in interfero- 
metric observations. We focus our study on the effects produced 



by turbulence in the dynamic range and astrometric accuracy 
after a phase-referenced calibration between a strong (calibra- 
tor) source and a weak source, located a few degrees away. This 
study is numerically equivalent to the study of the deformation 
of a wide-field interferometric image at any point located at a 
given distance from the center of the field (i.e., the phase cen- 
ter of the image). In both cases, the phases introduced by the 
atmosphere in the signal of each antenna for the different point- 
ing directions are the same, so the effects of the atmosphere in 
Fourier space (and therefore on the sky plane) will also be the 
same. 

The results here reported are an extended version of those 
previously reported in the SKA memo by Marti- Vidal et al. 
(2009). In the next section, we describe the details of the ar- 
ray distribution used, as well as the characteristics of the simu- 
lated observations. In Sect. [3] we describe how the noise from 
the atmosphere and the receivers was added to the visibilities 
and in Sect. |4] describe the procedures followed in our Monte 
Carlo analysis. In Sect. [5] we present the main results obtained; 
in Sect. [6] we summarize our conclusions. 



2. Array geometry and sensitivity 

We simulated an interferometric array similar to the planned sta- 
tion distribution of the Square Kilometre Array (SKA). We sim- 
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ulated a total of 200 stations distributed in the following way: 
50% are randomly distributed within a circle of 5 km radius (in- 
ner core); 25% are distributed outside this circle up to a distance 
of 150km (core), following 5 equiangular spiral arms; the re- 
maining antennae are distributed following the same spiral arms, 
but up to a distance of 3000 km from the inner core. This ar- 
ray distribution is similar to that used in Vir Lai, Lobanov & 
Jimenez-Monferrer (2009). The curvature of the Earth surface 
was taken into account in our simulations. We show the result- 
ing array distribution in Fig. Q] We also repeated all the simu- 
lations here reported, but subtracting a subset of 100 (randomly 
selected) stations from the array, to check the sensitivity of the 
main conclusions of this paper on different array distributions 
(see Appendix lAl. We also show in Fig. Q] the modified array 
after subtracting the 100 stations. 

2.1. Sensitivity and bandwidth 

We simulated interferometric observations using 16 different 
frequencies, which span in logarithmic bins from 150 MHz to 
24 GHz (this is the theoretical frequency window of the SKA). 
According to Jones (2004), the maximum observing bandwidth 
of the SKA will be around 25% of the central observing fre- 
quency (up to a maximum bandwidth of 4 GHz for all fre- 
quencies above 16 GHz). This (maximum) frequency-dependent 
bandwidth translates in our simulations into a changing sensitiv- 
ity of the SKA as a function of frequency. 

The sensitivities of the simulated stations were also chosen 
to be similar to those of the SKA, which were taken from Jones 
(2004). These values are set for an elevation of 45 degrees and 
differ from those given in Schilizzi et al. (120071 ). but the use of 
the values given in Schilizzi et al. (2007 ), instead, does not affect 
the main conclusions of this paper. We interpolated the sensitiv- 
ities given in Table 1 of Jones (2004) to the frequencies used in 
our simulations. In Fig. [2] we show the station sensitivities used. 
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Fig. 2. Station sensitivities (i.e., effective areas over system tempera- 
tures) used in our simulations. 



2.2. Source position 

We set the target source coordinates at the zenith of the array 
center and the calibrator at an hour angle of degrees, also with 
respect to the array center. This position of the sources mini- 
mizes the optical paths of the signals through the atmosphere 
(since both, source and calibrator, are at maximum elevations), 
thus enhancing the quality of the phase-referenced observations. 



The results given in this paper should be interpreted according 
to this issue. 

If the source would be located far from the zenith, the map- 
ping function of the tropospheric delay and the finite width of the 
ionosphere would increase the effect of turbulence on the phase- 
referenced visibilities of the target. Additionally, the uv coverage 
of the interferometer would have shorter projected baselines in 
declination, thus decreasing the synthesized resolution in decli- 
nation. Therefore, we want to stress that setting the sources at 
maximum elevations is a key limiting factor of the simulations 
here reported, especially if they are to be compared to real ob- 
servations. 



3. Noise model 

We simulated phase-referencing observations in the following 
way: we assumed that the calibrator source is sufficiently strong 
to allow for a perfect antenna-gain calibration at its location; we 
then determined the image of the target source by computing the 
differential antenna-gain errors expected at the target location. 
Therefore, under the effect of atmospheric turbulence, these re- 
sults depend on the calibrator-to-target separation. 

We implemented two kinds of atmospheric turbulence. The 
first turbulence was associated to the ionosphere (the free elec- 
tron content, which introduces dispersion in the radiation) and 
the other turbulence was associated to the wet troposphere (the 
water vapour, close to the earth surface, which is in a state of 
no thermodynamic equilibrium). The effect of ionospheric tur- 
bulence on the signal phase varies as v _1 , affecting the low- 
frequency observations; the effect of the wet troposphere on the 
phase varies as v, affecting the high-frequency observations. The 
dry troposphere (which is more homogeneously distributed over 
each station than the wet troposphere) was not considered in our 
simulations, since it can be easily modeled and removed from 
the data to a level lower than the effects coming from the wa- 
ter vapor and the ionosphere. Models of the turbulence from the 
ionosphere and troposphere can be found in many publications 
(e.g., Thomson, Moran, & Swenson |1991| ). Here, it is suffice to 
say that this turbulence follows a Kolmogorov distribution. This 
distribution has a phase structure function given by 

D^e) = < (O(0 O ) - Wo + e)f > oc e 5 ' 3 (l) 

where O(#o) is the phase added by the turbulent screen to the sig- 
nal of a source located at 6q. The brackets < ... > represent aver- 
aging over all pointing directions located at a distance 8 from the 
point at Go. The Kolmogorov distribution is fractal-like, so both, 
ionosphere and wet troposphere, have essentially the same phase 
distribution, despite of a global scaling factor between them. 

The global factors for both distributions (ionosphere and 
troposphere) were computed according to the typical values of 
ionospheric and tropospheric conditions. For the ionosphere, the 
Fried length (i.e., distance in the ionosphere for which the struc- 
ture function rises to lrad 2 ) was set to 3 km at 100 MHz. For 
the wet troposphere, we set the parameter C 2 n L (i.e., the inte- 
gral of the profile of C\ along the zenith direction) to 10" m 1 ' 3 
(vid. Eq. 13.100 and Table 13.2 of Thomson, Moran, & Swenson 
11991k this value translates into a Fried length of 3 km for a fre- 
quency of ~ 22 GHz. Since the Kolmogorov distribution is self- 
similar, it is possible to adapt the results here reported to any 
other atmospheric conditions (see Sect. 15.3b . just accordingly 
scaling the source separation to the Fried length of the iono- 
sphere (for low-frequency observations) or the wet troposphere 
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Fig. 1. Array distribution used in our simulations (empty and filled squares). Axes are relative Longitude (horizontal axis) and 
Latitude (vertical axis) in km. Left, the whole array. Center, a zoom to the core. Right, a zoom to the inner core. The stations marked 
with empty squares were removed from the array in a second run of our simulations, to check the dependence of the results on 
different array distributions (see text). 



(for high-frequency observations). We must notice that the self- 
similarity of the tropospheric turbulence does not hold for very 
large scales (the typical baseline lengths in VLBI observations), 
since there is a saturation in the power spectrum of the distribu- 
tion (see, e.g., Thomson, Moran, & Swenson 1991). However, 
this is not important in our analysis, since we did not use the 
absolute phase of the signal coming from a given direction in 
the sky, but computed the differential effects at each station from 
two different (closeby) directions, which depend on short-scale 
turbulence. Therefore, the saturation of tropospheric turbulence 
at large scales does not affect our results. 

We computed the differential effects from the turbulent at- 
mosphere in two ways. For the antennas of the core (within 
the central 300 km) we generated synthetic phase screens for 
the ionosphere and troposphere. We show an example of one 
such screen in Fig. [3] We notice that this figure could represent 
either ionospheric or tropospheric turbulence in our modeling, 
just by scaling the screen by the corresponding factor. Two dif- 
ferent screens were generated in each Monte Carlo simulation. 
The screen for simulating the ionosphere was put at a height of 
300 km and the screen for simulating the troposphere was put 
at a height of 5 km. For the antennas out of the core, we com- 
puted the term <S>(0q) - <E>(0o + S) separately. We proceeded this 
way (i.e., we generated a phase screen only for the core anten- 
nas, thus without generating a much larger screen for the whole 
array), because the distances between stations out of the core are 
large enough to ensure that the cross-correlation of turbulence 
above different stations is negligible compared to the correlation 
between those on the calibrator and target source for the same 
station. This numerical strategy also speeded up our simulations. 

It must be noticed that we did not introduce any time evolu- 
tion of the turbulent phase screens in our simulations. Any evolu- 
tion of the turbulence could dramatically affect the observations 
if the acquisition times were larger than the coherence time of the 
signal, which depends on the evolution of the turbulence and the 
observing frequency. However, for snapshot-like observations, 
of the order of a fraction of a minute or so, we could consider, as 
a good first approximation, a constant turbulence phase screen. 

Noise from the receivers was added to our model by gener- 
ating a random Gaussian noise in the real and imaginary parts 
of the visibilities. The mean deviation, <x, of the Gaussian noise 
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Fig. 3. Example of a turbulent phase screen with Kolmogorov 
statistics. The grey scale shows variations of optical-path phases, 
normalized between -1 and 1 radians. The final values of the 
phases depend (given the self-similarity of the distribution) on 
a global factor related to the observing frequency and the iono- 
spheric and/or wet tropospheric conditions. 



added to the visibilities was (e.g., Thomson, Moran, & Swenson 
fTWTI Eq. 6.43): 



i] Q VAvAr Sa 
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where k is the Boltzmann constant, tjq is the relative loss of sig- 
nal due to the correlator quantization (we used jjq = 0.5), Av is 
the observing bandwidth, Af is the observing time, and Sa is the 
sensitivity of the stations (collecting area over system tempera- 
ture, shown in Fig.[2]l. 



4. Estimate of dynamic range and astrometric 
precision 

We simulated different sets of phase-referenced observations. 
In all cases, the observations were snapshots with a duration 
to - 60 s. Longer observing times, t, would, in principle, in- 
crease the dynamic ranges and astrometric precisions shown in 
all the following sections as yftftQ, as long as the changing atmo- 
sphere (and, therefore, the changing source positions and shapes) 
would not introduce important smearing effects in the images af- 
ter the combination of all visibilities. 

In a first run of simulations, we generated visibilities of tar- 
gets with flux densities of 0.1, 1, and 10/iJy with a separation of 
5 degrees between target and calibrator. A total of 1500 sim- 
ulations were performed for each flux density and frequency. 
We used such a large separation between calibrator and target, 
because these simulations of phase-referenced observations can 
also be applied to the study of deformations of wide-field images 
under the effects of a turbulent atmosphere. 

In a second run of simulations, we studied the effects of 
the atmosphere as a function of calibrator-to-target separation. 
For that purpouse, we simulated 1500 observations at 1420 MHz 
(i.e., the Hydrogen line) of a source with 1 //Jy for different sep- 
arations from the calibrator (2, 3, 4, 5, and 6 degrees). 

In a third run of simulations, we used only one Kolmogorov 
screen (which can represent either ionospheric or tropospheric 
turbulence, depending on the observing frequency) with differ- 
ent Fried lengths, to study the scalability of the simulations for 
different source separations and/or atmospheric conditions. 

In all these simulations, we added the noise from the at- 
mosphere and the noise from the receivers. For each simulated 
phase-reference image, obtained by applying uniform weighting 
to the visibilities, the brightness peak was found and the corre- 
sponding point source was subtrated from the visibilities. For the 
substraction of the point source, the brightness peak was shifted 
to the phase center of the image by multiplying the visibilities 
by the corresponding plane-wave factor in Fourier space. Then, 
the flux density of the point source was estimated as the aver- 
age of the real part of the resulting visibilities, and the resulting 
point-source model was substracted from the data. Afterwards, a 
Fourier inversion of the new visibilities resulted in the image of 
residuals, from which the root-mean-square (rms) of all the pix- 
els was computed. On the one hand, the deviation of the bright- 
ness peak with respect to the image center was taken as the as- 
trometry error of that image. On the other hand, the source peak 
divided by the rms of the residuals was taken as the dynamic 
range. In Fig. [4] we show the distribution of astrometric devia- 
tions and dynamic ranges for the case of a target source of 1 /iJy 
observed at 1420 MHz (which corresponds to an interferometric 
beam of ~ 13 mas) located at 5 degrees from the calibrator. Once 
the distributions like those shown in Fig. [4] were obtained, we 
computed the standard deviation of astrometric corrections and 
the mean value of dynamic ranges for each source flux density, 
frequency, and separation. The first quantity was our estimate 
of the astrometric uncertainty, and the second quantity was an 
estimate of the achievable dynamic range. 
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Fig. 4. Distribution of right ascension shifts (a) and dynamic 
ranges (b) of 1500 snapshot phase-referenced images, simulated 
at 1420 MHz, for a 1 yt/Jy target source located at 5 deg from its 
calibrator. The calibrator source is located at an hour angle of 
and the target source is located at the zenith of the array center. 



5. Results 

5. 1 . Observing frequency and signal decoherence 

If the atmospheric turbulence is not taken into account and only 
the noise from the receivers is added to the visibilities, our sim- 
ulations reproduce the dynamic ranges given by Eq. 6.53 of 
Thomson, Moran, & Swenson (1 1 99 1 I t. as expected. Additionally, 
the noise from the receivers does not introduce considerable 
changes in the source position of the phase-referenced images 
(changes of the order of lOyuas or lower). 

When the turbulent ionosphere and wet troposphere are 
added to the simulations, the dynamic range of the images is no- 
tably affected, especially at low (< 1 GHz) and high (> 10 GHz) 
frequencies. In Fig. [5] we show phase-referenced images of a 
1 fiJy source, located at 5 deg. from its calibrator, observed at 
0.5, 5, and 15 GHz. It can be readily seen that the addition of 
effects coming from the atmospheric turbulence maps into an 
important extra noise in the images at 0.5 and 15 GHz, but not at 
5 GHz. 

Following the algorithm described in Sect. [4] we obtained 
the astrometric uncertainties and dynamic ranges shown in Fig. 
[6] For very low frequencies (below ~500 MHz) the ionosphere 
prevents a clear and precise detection of all sources, no matter 
their flux densities. For higher frequencies, the astrometric un- 
certainty decreases notably (mainly because of the dependence 
of ionospheric effects with v -1 ) and gets limited only by diffrac- 
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Fig. 5. Simulated phase-reference images of a 1 //Jy source located at 5 deg from its calibrator. Noise from the receiving system 
(top) and, additionally, noise due to the turbulent atmosphere (bottom) have been added to the visibilities. 



tion and sensitivity between 1 and 10 GHz (this frequency win- 
dow slightly depends on the source flux density, as it can be seen 
in the figure). For higher frequencies, the wet troposphere begins 
to affect the astrometric uncertainty, which rises up to around 
10 mas for the highest frequencies. We find that the best astro- 
metric accuracy, at least for reasonably well-detected sources, is 
achieved for frequencies around 4 GHz. This is where the iono- 
spheric and (wet) tropospheric components are roughly equal. 



The dynamic range of the phase-referenced images is highly 
limited by the atmosphere. When the atmosphere adds noise to 
the visibility phases, there is an extra rms added to the residual 
images, which depends on the visibility amplitudes, thus limit- 
ing the achievable dynamic range no matter the flux density of 
the source; that is, if the source flux density is higher, the noise 
of the image will also be higher. This limitation is, of course, 
more important for the brightest sources. In our case, the bright- 
est source has a flux density of 10 /Jy. For this source, the max- 
imum dynamic range achieved is only 110, which is ~ 30 times 
smaller than the dynamic range that would be obtained without 
the atmosphere. This situation can be also understood in another 
way: the rms of the final image is divided into two components, 
which are added in quadrature. One component, cr,/,, comes from 
the receiver noise and is independent of the source flux den- 
sity. The other component, cr at , comes from the atmospheric re- 
fraction and is equal to a percentage of the source flux density 



(cr at — K r S , where S is the source flux density and K r depends 
on the atmospheric refraction). Hence, the dynamic range, D, is 



D 



4 



crt, + cr„ 



4 



(3) 



For large flux densities (S » cr,/,), the achievable dynamic 
range will saturate to a value dependent on the atmospheric con- 
ditions (i.e., D — ♦ 1/K r ) and independent on the source flux den- 
sity and the sensitivity of the stations. As seen in Fig. [6] despite 
of the large difference between the simulated flux densities (1 
and 10 //Jy), the dynamic range saturates to a value of ~100. 

5.2. Angular separation and signal decoherence 

The results shown in the previous subsection correspond to a 
separation of 5 degrees between source and calibrator. These re- 
sults change when the angular separation changes. We computed 
astrometric uncertainties and dynamic ranges for a source with 
a flux densitity of 1 fiJy located at 2, 3, 4, 5, and 6 degrees from 
its calibrator. Noise from the atmosphere and the receivers was 
taken into account in these simulations. We used an observing 
frequency of 1420 MHz (the Hydrogen line) which is inside the 
frequency window where the atmospheric effects are minimised. 
Therefore, all the astrometric errors derived were small (of the 
order of a few mas), allowing us to use image sizes small enough 
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Fig. 6. Astrometric accuracy (a) and dynamic range (b) as a func- 
tion of frequency, for the case with atmospheric turbulence and 
a separation of 5 degrees between calibrator and target. The cali- 
brator source is located at an hour angle of and the target source 
is located at the zenith of the array center. Different lines corre- 
spond to different target source flux densities (10 fi Jy, continuous 
line; 1 /iJy, dashed line; 0.1 fiJy, dotted line). 
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Fig. 7. Dots: simulated astrometric accuracy (a) and dynamic 
range (b) as a function of angular separation between calibra- 
tor and source, for observations at 1420 MHz and a target flux 
density of 1 /iJy. Lines: proposed phenomenological models. 



to sample the beam with more pixels (~30 pixels) using a grid 
of 1024x1024 pixels. This fine gridding of the beam allowed for 
a more accurate determination of the location of the image peak 
and, therefore, a better estimate of the astrometric error. The re- 
sults obtained are shown in Fig. [7] In that figure, we also plot two 
analytical (phenomenological) models for the estimate of the in- 
crease of astrometric uncertainty and the loss of dynamic range 
(i.e., degree of signal decoherence) as a function of angular sep- 
aration. On the one hand, the phenomenological model proposed 
for the estimate of loss of dynamic range is: 



D 



Do 



(4) 



where D is the dynamic range, 8 is the angular separation be- 
tween target and calibrator, Do is the dynamic range without at- 
mosphere (i.e., when the calibrator-to-target separation, 6, tends 
to 0), and k and (3 are two parameters related to the atmospheric 
conditions, source flux density, and observing frequency. This 
equation is just Eq. [3] but setting K r = k O 8 ^ 2 . 

As it can be seen, this model fits well to the simulations. We 
obtain k = (2.5 + 0.2) x 10~ 5 deg^ and/3 = 1.51 ±0.06. 

On the other hand, the proposed phenomenological model 
for the increase of astrometric uncertainty is: 



cr = — - (1 + k 1 1 
D 



(5) 



where <x is the astrometric uncertainty, <x is the diffraction limit 
(i.e., half the size of the beam), and k' and /?' are two parame- 
ters also related to the atmospheric conditions, source flux den- 
sity, and observing frequency. We fit k' = 0.57 ± 0.04 deg^ and 
/3' = 1.02 ± 0.05. From a different approach in the treatment of 
the noise from the atmosphere and limited to VLBA and EVN 
arrays, Pradel, Chariot, & Lestrade (2006) obtained a similar de- 
pendence of <j with 9. 

We notice that if we change Do by D in Eq. [5] the new fit- 
ted k' and 0' are 1.03+0.09 and 0.07+0.06, respectively. This 
new value of /3' is compatible with zero. In other words, the 
diffraction limit divided by the dynamic range of the image is 
an excellent estimator of the astrometric uncertainty, at least for 
the range of simulated calibrator-to-target separations at 1 .4 GHz 
(which falls within the frequency window where the atmospheric 
effects are minimized). 

For calibrator-to-target separations larger than ~6 degrees, 
the situation changes. We simulated phase-referenced images for 
calibrator-to-target separations up to 12 degrees, and found that 
the model of dynamic range given by Eq.|4]is still valid, but the 
astrometric uncertainty increases faster, with /?' = 2.38 + 0.08 
(J3' = 1.38, if we change D by D in Eq.g). This last /?' fits well 
to the astrometric uncertainties for large source separations, but 
the fit is worse for separations smaller than 5-6 degrees. 



I. Martf-Vidal et al.: Atmospheric turbulence in interferometric images 



7 



5.3. Scalability of the results and use of turbulence models in 
the data calibration 

In the previous subsections, we report on the effects of atmo- 
spheric turbulence in phase-referenced (and wide-field) inter- 
ferometric images using fixed values for the Fried lengths of 
the Kolmogorov distributions of the ionosphere and wet tropo- 
sphere. Since the Kolmogorov distribution is self-similar, the re- 
sults reported can be scaled and adapted to other atmospheric 
conditions. Indeed, these simulations can also be used to esti- 
mate the limiting dynamic range and astrometric uncertainty if 
an a priori model of the tropospheric and/or ionospheric turbu- 
lence is used in the imaging. In these cases, the effective Fried 
length, r e f, to compare to our simulations can be estimated as 



a wide-angle coverage at very low (or high) frequencies. For in- 
stance, an image of 10x10 degrees with a dynamic range of ~10 6 
at a frequency of 500 MHz (which translates into rj ~ 5 x 10~ 5 ) 
would require, for the ionospheric screen used in the previous 
subsections, a model of the ionospheric turbulence distribution 
with an uncertainty lower than 3 x 10~ 5 (observing time of 1 
minute). If the observing time increases to 10 hours, the mini- 
mum required uncertainty in the model of the ionosphere would 
increase to 7 x 10~ 4 , but such a high precision should be kept 
during the whole set of observations. If the dynamic range de- 
creases to 10 5 , the required minimum uncertainty for the iono- 
sphere model would decrease to 7 x 10~ 3 (i.e., 0.7%) for an ob- 
serving time of 10 hours. 



r ef 



\<t> ~ <Pmod\ 



> 



(6) 



where ro is the Fried length of the real turbulence and the other 
factor is related to the fractional precision of the turbulence 
model: (p mo d is the phase computed from the turbulence model 
at a given point in the sky and (f> is that corresponding to the real 
turbulence; the brackets < ... > represent averaging over the field 
of view. If the a priori model of the ionospheric electron distri- 
bution is accurate to a given precision level, the effective Fried 
length to use for the ionosphere will be that corresponding to the 
residual turbulence (i.e., the difference between the model and 
the real turbulence). For instance, in the case of a model of the 
ionospheric electron content with a 99% accuracy, the effective 
Fried length will be r e f = 100 ro; if the accuracy increases to 
99.9%, r e f = 1000 ro. In Fig. [8] we show D max , the maximum 
dynamic range (i.e., for a source with an infinite flux density, so 
cr,h = in Eq. as a function of parameter rj, which we define 
as 
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Fig. 8. Maximum dynamic range as a function of 77 (see Eq. |7]i 
and for a 1-scan snapshot of 60 s. 



h sinfl 



77 = 



(7) 



In this equation, h is the height of the phase screen and 9 is 
the calibrator-to-target separation (or half the size of the wide- 
field image). Equation Q can be used, together with Fig. [8] to 
compute the maximum achievable dynamic range for many dif- 
ferent combinations of source separations, atmospheric condi- 
tions, and observing frequencies {r e f oc y for the ionosphere and 
r e f oc y _1 for the trosposphere). We notice, however, that Fig. 
[8] has been generated using only one Kolmogorov screen, so it 
is applicable to ionospheric dispersion (for low frequencies) or 
tropospheric dispersion (for high frequencies), but not to a situ- 
ation where ionospheric and tropospheric effects are similar. In 
these cases, and as a first approximation, we could set 
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Fig. 9. Peak flux density of the phase-referenced image, relative 
to the real peak flux density of the source, as a function of 77 (see 
Eq.|7]i and for a 1-scan snapshot of 60s. 



where r lo „ and r trop are the (effective) Fried lengths of the iono- 
sphere and troposphere, respectively, and h wn and h rrop are the 
heights of each phase screen. Also shown in Fig.[8]is the fitting 
model 

D max = D\ rf' (9) 

where £>i = 48.38 ± 0.10 and /3i = -1.000 ± 0.002. Figure[8] 
indicates how difficult is to obtain a high-contrast image with 



For a more comprehensive representation of our results, we 
show in Fig. [10] how the achievable dynamic range (computed 
from Eq. [9} depends on the uncertainty in the model of atmo- 
spheric turbulence used in the data calibration. We show this 
relationship for different observing frequencies and calibrator- 
to-target separations. For instance, a dynamic range 10 4 in ob- 
servations at 100 MHz for a calibrator-to-target separation of 1 
deg. (i.e., the same for a wide-field image of 2 x 2 deg.) would 
require a turbulence model with an accuracy of ~ 99.9% for an 
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observing time of 60s. This requirement would soften to an ac- 
curacy of 97 - 98% for an observing time of 6000 s (provided 
the dynamic range increases as the square root of the observing 
time). 
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Fig. 10. Maximum dynamic range as a function of the uncer- 
tainty in the model of the atmospheric turbulence for a 1-scan 
snapshot of 60 s. 

For completeness, we also computed the loss of recovered 
flux density of the source due to the turbulent atmosphere. In Fig. 
[9] we show the ratio of peak flux densities between the phase- 
referenced images and those computed without the effects from 
atmospheric turbulence. For dynamic ranges of 40-50, the loss 
of flux density can be as large as 25%. The phenomenological 
model (also shown in the figure) used to fit the data is 



where D 2 = 0.32 ± 0.02 and/? 2 = 1.87 + 0.01. 

Applicability of Eq. [9] and [10] is not restricted to the array 
used in the simulations here reported. The exponents j3\ and 
$2 only depend on the structure of the atmospheric turbulence 
and are thus independent of the interferometer used in the ob- 
servations. However, the parameters D\ and D2 also depend on 
the stations of the interferometer. Therefore, Eq. [9] and [10] can 
be adapted to any other interferometer by finding the right val- 
ues of D\ and Di. As an example of this generalization of Eq. 
[9] and [10] Martf- Vidal et al. (120101) have studied the achiev- 
able dynamic range in phase-referenced observations with the 
Very Long Baseline Array (VLBA) at 8.4 GHz and 15 GHz. This 
study has been performed from quasi-simultaneous observations 
of 13 sources located at separations ranging from 1.5 to 20.5 
degrees. These authors have been able to model the dynamic 
ranges obtained in the phase-referenced images and the loss of 
recovered flux densities (phase-referenced images compared to 
images obtained from self-calibrated visibilities) using Eq.[9]and 
[10]with values for D\ and D2 different to those here reported, but 
using the same values here reported for the exponents f3\ and y02- 

6. Conclusions 

We report on Monte Carlo estimates of the sensitivity and astro- 
metric precision of an interferometric array, with a station dis- 
tribution similar to that of the planned SKA, as a function of 
observing frequency, flux density, and source separation. These 
results can also be applied to other array distributions by taking 



into account the corresponding correction factors. Our estimates 
are based on simulations of snapshot phase-referenced observa- 
tions, in which we take into account several effects from the tur- 
bulent atmosphere and the finite temperature of the receivers. 
We find that the astrometric uncertainty strongly depends on the 
observing frequency and smoothly increases as the source sepa- 
ration increases. For frequencies below ~1 GHz, ionospheric ef- 
fects dominate and the astrometry uncertainties (when the source 
is detectable) can be as large as ~1 as. For frequencies between 1 
and 10 GHz (these values slightly depend on the source flux den- 
sity) atmospheric effects are minimum and we roughly reach the 
theoretical precision of the interferometer. Above these frequen- 
cies, the wet troposphere begins to dominate and the astrometric 
uncertainty increases to ~ 10 mas for the highest simulated fre- 
quency (25 GHz). The dynamic range of the images is strongly 
limited by atmospheric turbulence at all frequencies and for all 
flux densities (it can decrease, in the worse cases, several orders 
of magnitude). 

We propose analytical models for the loose of dynamic range 
and astrometric accuracy as a function of distance between cal- 
ibrator and target source. These expressions could also be used 
to estimate the deformations and local dynamic ranges of wide- 
field images as a function of distance to the image phase center 
(i.e., the point in the sky where the data correlation is centered). 
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Appendix A: Complementary simulations: different 
number of stations and array sensitivities 

Our simulations are based on a given realization of the SKA. 
However, the main structure of the array distribution used in our 
simulations is not exclusive of the SKA. Other interferometric 
arrays, like ALMA or LOFAR, are being built with similar sta- 
tion distributions, consisting on a compact core and several ex- 
tensions with the shape of spiral arms. Hence, our study can be 
extended to those arrays by taking into account the difference 
between the number of stations and the station sensitivities. 

The amount of noise added to the data is proportional to the 
number of stations, since each station receives the signal through 
the turbulent atmosphere. However, for cases of clear source de- 
tections (D > 20-30), the dynamic range does not depend (or the 
dependence is weak) on the thermal noise of the receivers (see 
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Sect. 15. U . Therefore, to estimate the achievable dynamic range 
for an array with a different number of stations, the results shown 
in Fig. [6] should be divided by N/N', where N is the number of 
stations used in our simulations (N = 200) and N' is the number 
of stations of the other array. This is true for detections with a 
relatively large dynamic range (D > 20 - 30). For weak sources, 
the noise from the receivers may also contribute to the rms of the 
residual images, so the factor to apply in these cases should be 
(Ncr')l(N'cr), where <x is the thermal noise of the stations used 
in our simulations and cr' is that of the other array. 

We repeated the simulations described in Sect. l5.1l using dif- 
ferent arrays to compare the results with those obtained with the 
original array. On the one hand, we created a smaller array by 
subtracting 100 stations (those marked with empty squares in 
Fig. [1} from the original array. On the other hand, we created 
another array with all the 200 stations, but decreasing their sen- 
sitivity by a factor 2. We show the results obtained in Fig. IA. II 
Special care must be taken in the interpretation of these figures, 
given that the computed ratios of dynamic ranges are only mean- 
ingful when the detections of the sources are clear (i.e. when 
no spurious noise peaks appear stronger than the source). This 
is true for D>20-30, which approximately corresponds to fre- 
quencies between 1 and 10 GHz (although it slightly depends on 
the source flux density, see Fig. [6j>. 

With these considerations into account, we find that the ra- 
tio of dynamic ranges for the array with 100 stations falls be- 
tween 0.7 and 0.5 compared to the array with 200 stations. The 
expected value is 0.5 (since N' = 0.5N and cr' = cr). Other fac- 
tors, like the different coverages of Fourier space by both arrays, 
could be affecting the dynamic range of the images, thus increas- 
ing the ratio in some cases. For the array with lower station sen- 
sitivities (but the same number of stations), the ratio of dynamic 
ranges falls between 0.75 and 1 for the strongest sources (as ex- 
pected, since N' - N and the noise from the receivers is much 
smaller than the noise from the atmosphere), while are close to 
0.5 for the weakest source (also as expected, since the thermal 
noise from the receivers begins to dominate in this case). 
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Fig. A.l. Ratios of dynamic ranges obtained with our original ar- 
ray and those obtained with a subarray of 100 stations (a) and an 
array with half the sensitivity of the original array (b). Different 
lines correspond to different target source flux densities (lOjuJy, 
continuous line; 1 ptJy dashed line; 0.1 ptJy, dotted line). 



